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Before corresponding points in images taken with two cameras can be 

distances to objects in a scene, one has to determine the position and 

of one camera relative to the other. This is the classic photogrammetric 

relative orientation, central to the interpretation of binocular stereo 

Iterative methods for determining relative orientation were developed 

without them we would not have most of the topographic maps we do 

ative orientation is also of importance in the recovery of motion and 

an image sequence when successive frames are widely separated in time. 

motion vision are rediscovering some of the methods of photogrammetry. 

here is a particularly simple iterative scheme for recovering relative 

that, unlike existing methods, does not require a good initial guess for 

and the rotation. The data required is a set of pairs of corresponding 

two projection centers to points in the scene. It is well known that at 

of rays are needed. Less appears to be known about the existence of 

utions and their interpretation. These issues are discussed here in detail. 

ious determination of all of the parameters of relative orientation is 

when the observed points lie on a critical surface. These surfaces and 

erate forms are analysed here as well. 
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Relative Orientation 1 

1. Introduction 

The positions of corresponding points in two images can be used to determine 
the positions of points in the environment, provided that the position and 
orientation of one camera with respect to the other is known. Given the 
internal geometry of the camera, including its focal length and the location of 
the principal point, rays can be constructed by connecting the points in the 
images to their corresponding projection centers. These rays, when extended, 
intersect at the point in the scene that gave rise to the image points. This is 
how binocular stereo data is used to determine the positions of points in the 
environment after the correspondence problem has been solved. 

It is also the method used in motion vision when feature points are tracked 
and the image displacements that occur in the time between two successive 
frames are relatively large (see for example [Ullman 1979] and [Tsai & Huang 
1984]). The connection between these two problems has not attracted much 
attention before, nor has the relationship of motion vision to some aspects 
of photogrammetry (but see [Longuet-Higgins 1981]). It turns out, for exam- 
ple, that the well known motion field equations [Longuet-Higgins & Prazdny 
1980, Bruss & Horn 1983] are just the parallax equations of photogrammetry 
[Hallert 1960, Moffit & Mikhail 1980] that occur in the incremental adjust- 
ment of relative orientation. Most papers on relative orientation only give the 
equation for y-parallax, corresponding to the equation for the j/-component of 
the motion field (see for example the first equation in [Gill 1964], equation (1) 
in [Jochmann 1965], and equation (6) in [Oswal 1967]). Some papers actually 
give equations for both x- and y-parallax (see for example equation (9) in 
[Bender 1967]). 

In both binocular stereo and large displacement motion vision analysis, 
it is necessary to first determine the relative orientation of one camera with 
respect to the other. The relative orientation can be found if a sufficiently large 
set of pairs of corresponding image points have been identified [Thompson 
1959b, Thompson 1968, Ghosh 1972, Schwidefsky 1973, Slama et al. 1980, 
Moffit & Mikhail 1980, Wolf 1983, Horn 1986]. 

Let us use the terms left and right to distinguish the two cameras (in the 
case of the application to motion vision these will be the camera positions and 
orientations corresponding to the earlier and the later frames respectively). 1 
The ray from the center of projection of the left camera to the center of 
projection of the right camera is called the baseline (see Fig. 1). A coordinate 
system can be erected at each projection center, with one axis along the 

x In what follows we use the coordinate system of the right (or later) camera as 
the reference. One can simply interchange left and right if it happens to be more 
convenient to use the coordinate system of the left (or earlier) camera. 



optical axis, that is, perpendicular to the image plane. The other two axes 
are in each case parallel to two convenient orthogonal directions in the image 
plane (such as the edges of the image, or lines connecting pairs of fiducial 
marks). The rotation of the left camera coordinate system with respect to 
the right is called the orientation. 
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Figure 1. Points in the environment are viewed from two camera positions. 
The relative orientation is the direction of the baseline b, and the rotation 
R relating the left and right coordinate systems. The directions of rays 
to at least five scene points must be known in both camera coordinate 

systems. 

Note that we cannot determine the length of the baseline without knowledge 
about the length of a line in the scene, since the ray directions are unchanged 
if we scale all of the distances in the scene and the baseline by the same 
positive scale-factor. This means that we should treat the baseline as a unit 
vector, and that there are really only five unknowns — three for the rotation 
and two for the direction of the baseline. 2 

2. Existing Solution Methods 

Various empirical procedures have been devised for determining the relative 
orientation in an analog fashion. Most commonly used are stereoplotters, op- 



If we treat the baseline as a unit vector, its actual length becomes the unit of 
length for all other quantities. 
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tical devices that permit viewing of image pairs and superimposed synthetic 
features called floating marks. Differences in ray direction parallel to the 
baseline are called horizontal disparities (or ar-parallaxes), while differences in 
ray direction orthogonal to the baseline are called vertical disparities (or y- 
parallaxes). 3 Horizontal disparities encode distances to points on the surface 
and are the quantities sought after in measurement of the underlying topogra- 
phy. There should be no vertical disparities when the device is adjusted to the 
correct relative orientation, since the rays from the left and right projection 
center must lie in a plane that contains the baseline if they are to intersect. 

The methods used in practice to determine the correct relative orien- 
tation depend on successive adjustments to eliminate the vertical disparity 
at each of five or six carefully chosen points [Sailor 1960, Thompson 1964, 
Slama et al. 1980, Moffit & Mikhail 1980, Wolf 1974]. In each of the ad- 
justments a single parameter of the relative orientation is varied in order to 
remove the vertical disparity at one of the points. Which adjustment is made 
to eliminate the vertical disparity at a particular point depends on the partic- 
ular method is chosen. In each case, however, one of the adjustments, rather 
than being guided visually, is made by an amount that is calculated, using 
the measured values of earlier adjustments. The calculation is based on the 
assumptions that the surface being viewed can be approximated by a plane, 
that the baseline is roughly parallel to this plane, and that the optical axes of 
the two cameras are roughly perpendicular to this plane. The whole process 
is iterative in nature, since the reduction of vertical disparity at one point 
by means of an adjustment of a single parameter of the relative orientation 
disturbs the vertical disparity at the other points. Convergence is usually 
rapid if a good initial guess is available. It can be slow, however, when the 
assumptions on which the calculation is based are violated, such as in "acci- 
dented" or hilly terrain [Van Der Weele 1959-60]. These methods typically 
use Euler angles to represent three-dimensional rotations [Korn & Korn 1968] 
(traditionally denoted by the greek letters k, <f>, and a;). Euler angles have 
a number of short-comings for describing rotations that become particularly 
noticable when these angles become large. 

There also exist related digital procedures that converge rapidly when a 
good initial guess of the relative orientation is available, as is usually the case 
when one is interpreting aerial photography [Slama et al. 1980]. Not all of 
these methods use Euler angles. Thompson [1959b], for example, uses twice 
the Gibb's vector [Korn & Korn 1968] to represent rotations. These proce- 

3 This naming convention stems from the observation that, roughly speaking, in 
the usual viewing arrangement, horizontal disparities correspond to left-right 
displacements in the image, whereas vertical disparities correspond to up-down 
displacements. 



dures may fail to converge to the correct solution when the initial guess is far 
off the mark. In the application to motion vision, approximate translational 
and rotational components of the motion are often not known initially, so a 
procedure that depends on good initial guesses is not particularly useful. 

Normally, the directions of the rays are obtained from points in images 
generated by projection onto a planar surface. In this case the directions are 
confined to the field of view as determined by the active area of the image plane 
and the distance to the center of projection. The field of view is always less 
than a hemi-sphere, since only points in front of the camera can be imaged. 4 
The method described here applies, however, no matter how the directions 
to points in the scene are determined. There is no restriction on the possible 
directions of the rays. We do assume, however, that we can tell which of two 
opposite semi-infinite line-segments the point lies on. If a point lies on the 
correct line-segment we will say that it lies in front of the camera, otherwise 
it will be considered to be behind the camera. 

The problem of relative orientation is generally considered solved, and 
so has received little attention in the photogrammetric literature in recent 
times [Van Der Weele 1959-60]. In the annual index of Photogrammetric 
Engineering, for example, there is only one reference to the subject in the last 
ten years [Ghilani 1983] and six in the decade before that. This is very little 
in comparison to the large number of papers on this subject in the fifties, 
as well as the sixties, including [Gill 1964], [Sailor 1965], [Jochmann 1965], 
[Ghosh 1966], [Forrest 1966] and [Oswal 1967]. 

In this paper we discuss the relationship of relative orientation to the 
problem of motion vision in the situation where the motion between the ex- 
posure of successive frames is relatively large. Also, a new iterative algorithm 
is described here, as well as a way of dealing with the situation when there 
is no initial guess available for the rotation or the direction of the baseline. 
The advantages of the unit quaternion notation for representing rotations are 
illustrated as well. Finally, we discuss critical surfaces, surface shapes that 
lead to difficulties in establishing a unique relative orientation. 

3. Coplanarity Condition 

If the ray from the left camera and the corresponding ray from the right 
camera are to intersect, they must to lie in a plane that also contains the 
baseline. Thus, if b is the vector respresenting the baseline, r r the ray from 
the right projection center to the point in the scene and r/ the ray from the 



The field of view is larger than a hemi-sphere in some fish-eye lenses, where there 
is significant radial distortion. 
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left projection center to the point in the scene, then the triple product 

[b r} r r ] 
equals zero, where rj = Rot(r;) is the left ray rotated into the right camera's 
coordinate system. 5 This is the coplanarity condition (see Fig. 2). 







Figure 2. Two rays approach closest where they are intersected by a 
line perpendicular to both. If there is no measurement error, and the 
relative orientation has been recovered correctly, then the two rays actually 
intersect. In this case the two rays and the baseline lie in a common plane. 



We obtain one such constraint from each pair of rays. There will be an infinite 
number of solutions for the baseline and the rotation when there are fewer 
than five pairs of rays, since there are five unknowns and each pair of rays 
yields only one constraint. Conversely, if there are more than five pairs of 
rays, the constraints are likely to be inconsistent as the result of small errors 
in the measurements. In this case, no exact solution of the set of constraint 
equations exists, and it makes sense instead to minimize the sum of squares of 
errors in the constraint equations. In practice, one should use more than five 
pairs of rays in order to reduce the influence of measurement errors [Jochmann 
1965, Ghosh 1966]. We shall see later that the added information also allows 
one to eliminate spurious solutions. 



'The baseline vector b is here also measured in the coordinate system of the right 
camera. 
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4. What is the Appropriate Error Term? 

The triple product [b rj r r ] is zero when the left and right ray are coplanar 
with the baseline. It is not immediately apparent, however, that the triple 
product itself is necessarily the ideal measure of departure from best fit. It is 
worthwhile exploring the geometry of the two rays more carefully. Consider 
the points on the rays where they approach each other the closest (see Fig. 2). 
The line connecting these points will be perpendicular to both rays, and hence 
parallel to (rj x r r ). As a consequence, we can write 

ar', + 7(r{ x r r ) = b + /?r r , 

where a and ft are proportional to the distances along the left and the right ray 
to the points where they approach most closely, while 7 is proportional to the 
shortest distance between the rays. We can find 7 by taking the dot-product 
of the equality above with rj x r r . We obtain 

7||r',xr r || 2 = [br; r r ]. 

Similarly, taking the dot-products with r r x (r{ x r r ) and rj x (rj x r r ), we 
obtain 

a ||r} x r r || 2 = (b x r r ) • (r{ x r r ), 

/9||r , I xr r || 2 = (bxrJ)-(pJxr P ). 
The magnitudes of a and ft are the distances along the rays to the point of 
closest approach when r r and r{ are unit vectors. It turns out, however, that 
we are more interested here in the signs than in the magnitudes of a and ft. 

Normally, the points where the two rays approach the closest will be in 
front of both cameras, that is, both a and ft will be positive. If the estimated 
baseline or rotation is in error, however, then it is possible for one or both of 
the calculated parameters a and ft to be negative. We will use this observation 
later to distinguish amongst different apparent solutions. 

We have shown that the perpendicular distance between the left and the 
right ray is equal to ratio of the triple product [b r{ r r ] to the magnitude 
squared of (r{ x r r ). But the measurement errors are in the image, not in 
the scene. Thus a least-squares procedure should be based on an error in 
determining the direction of the rays, not on the distance of closest approach. 
To arrive at such a measure, we can look at the angle, 8 say, between the 
projections of the left and right ray into a plane perpendicular to b (see 
Fig. 3). This angle will be zero when the vertical disparity is zero, that is, 
when the left and right rays are coplanar with the baseline. 

The projections of r{ and r r into a plane perpendicular to b are given by 

r} = rj - (r' ; • b)b = (b x rj) x b, 
r r = r r — (r r • b)b = (b x r r ) x b, 



Relative Orientation 




Figure 3. One measure of the departure from the coplanarity condition 
is the angle 6 between the planes formed by the left ray and the baseline 
and the plane formed by the right ray and the baseline. The angle may 
be found by projection the two rays onto a plane perpendicular to the 
baseline. 



where we have used the fact that b is a unit vector. It is easy to show that 

Ur',11 = ||b x ri|| and ||r r || = ||b x r r || , 

keeping in mind again that b is a unit vector, and that b x rj and b x r r are 
perpendicular to b. The cross-product of the two projected vectors rj and r r 
will be parallel to the baseline and have magnitude proportional to the sine 
of the angle between the projected vectors, and so, 

[b r} r r ] _ [b rj r r ] 



sin# = 



\r',\ 



b x r{|| ||b x r r 



where we have used the fact that [b rj r r ] — [b r' { r r ], something that can 
easily be verified. 
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We could use the sine of the angle between the projected vectors directly 
as a measure of departure from best fit. This is not as good an idea as it may 
appear at first sight, because of what happens when one of the rays becomes 
nearly parallel to the baseline. In this case the angle will vary rapidly with 
small changes in the direction of the ray. Correspondingly, one of the terms 
in the denominator in the expression for the sine of the angle becomes small. 
It is better to normalize the expression by multiplying by the lengths of the 
projected vectors. Then we obtain the area of the parallelogram formed by 
the projections of the rays into a plane perpendicular to the baseline, namely, 

||f{||||p f .||8intf = [bf{p r ] = [bpjr r ]. 
This discussion confirms that the triple product itself is a good measure of the 
departure from best fit. This is convenient, since it makes the least squares 
analysis reasonably straightforward. If one desires to use a different error 
measure, one can weight the terms in the sums to follow accordingly. 

5. Least Squares Solution for the Baseline 

If the rotation is known, it is easy to find the best fit baseline, as we show next. 
This is useful, despite the fact that we do not usually know the rotation. The 
reason is that the ability to find the best baseline, given a rotation, reduces 
the dimensionality of the search space from five to three. 

Let {ri t i} and {r rj ,}, for i = l...n, be corresponding sets of left and 
right rays. We wish to minimize 

£ = D br t^] 2 = E( b 'Ki x M 2 . 

i=l i=l 

subject to the condition ||b|| = 1, where rj i is the rotated left ray r/^, as 
before. If we let c 2 =rJjX r r; ;, we can rewrite the sum in the simpler form 

25 = £> • cO 2 = b r (f>cr) b, 

1=1 \j=l / 

where we have used the equivalence b • c,- = b r Cj, which depends on the 
interpretation of column vectors as 3 x 1 matrices. The term Cicf is a dyadic 
product, a 3 x 3 matrix obtained by multiplying a 3 x 1 matrix by a 1 X 3 
matrix. 

The error sum is a quadratic form involving the real symmetric matrix 

C = £c,cf. 

t=l 

The minimum of such a quadratic form is the smallest eigenvalue of the matrix 

C, attained when b is the corresponding unit eigenvector (see, for example, 
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the discussion of Rayleigh's quotient in [Korn & Korn 1968]). This can be 
verified by introducing a Lagrangian multiplier A and minimizing 

E' = b r Cb + A(l-b T b), 
subject to the condition b T b = 1. Differentiating with respect to b and 
setting the result equal to zero yields: 

Cb = Ab. 
The error corresponding to a particular solution of this equation is found by 
premultiplying by b T : 

E = b T Cb = \b T b = X. 
The three eigenvalues of the real symmetric matrix C are non-negative, and 
can be found in closed form by solving a cubic equation, while each of the 
corresponding eigenvectors has components that are the solution of three ho- 
mogeneous equations in three unknowns [Korn & Korn 1968]. If the data 
are relatively free of measurement error, then the smallest eigenvalue will 
be much smaller than the other two, and a reasonable approximation to the 
sought-after result can be obtained by solving for the eigenvector using the 
assumption that the smallest eigenvalue is actually zero. This way one need 
not even solve the cubic equation (see also [Horn & Weldon 1988]). 

If b is a unit eigenvector, so is — b. Changing the sense of the baseline 
does not change the magnitude of the error term [b r{ r r ]. It does, however, 
change the signs of a, j3 and 7. One can decide which sense of the baseline 
direction is appropriate by determining the signs of ati and /% for i = 1 . . . n. 
Ideally, they should all be positive. The solution for the optimal baseline is 
not unique unless there are at least two pairs of corresponding rays. The 
reason is that the eigenvector we are looking for is not uniquely determined 
if more than one of the eigenvalues is zero, and the matrix has rank less than 
two if it is the sum of fewer than two dyadic products of independent vectors. 
This is not a significant restriction, however, since we need at least five pairs 
of rays to solve for the rotation anyway. 

6. Iterative Improvement of Relative Orientation. 

If one ignores the orthonormality of the rotation matrix, a set of nine homoge- 
neous linear equations can be obtained by a transformation of the coplanarity 
conditions that was first described in [Thompson 1959b]. These equations 
can be solved when eight pairs of corresponding ray directions are known 
[Longuet-Higgins 1981]. This is not a least-squares method that can make 
use of redundant measurements, nor can it be applied when fewer than eight 
points are given. The method is also strongly affected by measurement errors 
and fails for certain configurations of points [Longuet-Higgins 1984]. 
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No true closed-form solution of the least-squares problem has been found 
for the general case, where both baseline and rotation are unknown. However, 
it is possible to determine how the overall error is affected by small changes 
in the baseline and small changes in the rotation. This allows one to make 
iterative adjustments to the baseline and the rotation that reduce the sum of 
squares of errors. 

We can represent a small change in the baseline by an infinitesimal quan- 
tity 8b. If this change is to leave the length of the baseline unaltered, then 

lib + *b|| 2 = libit 2 , 



or 



|b|| 2 +2b^b + ||^b|| 2 = ||b|| 2 



If we ignore quantities of second-order, we obtain 

8b ■ b = 0, 
that is, 8b must be perpendicular to b. 

A small change in the rotation can be represented by an infinitesimal 
rotation vector 8u. The direction of this vector is parallel to the axis of rota- 
tion, while its magnitude is the angle of rotation. This incremental rotation 
takes the rotated left ray, rj, into 

r'/ = r{ + (6u x rj). 

This follows from Rodrigues' formula for the rotation of a vector r: 

cos 8 r + sin 8 (w x r) + (1 — cos 8){u ■ r) u>, 

if we let 8 — \\Su\\, w = 8u>/ ||<ta||, and note that 8u is an infinitesimal 
quantity. Finally then, we see that [b r\ r r ] becomes 

[(b + 6b) (rj + tfwxrj) r r ], 

or, 

[b rj r r ] + [6b r| r r ] + [b (Su x r',) r r ], 
if we ignore the term [^b {6lj x r{) r r ], because it contains the product of 
two infinitesimal quantities. We can expand two of the triple products in the 
expression above and obtain 

[b r{ r r ] + (r'j x r r ) • <5b + (r r x b) • (6u x rj), 

or 

t + c-6b + d-6u}, 

for short, where 

t = [b rj r r ], c = rj x r r , and d = rj x (r r x b). 
Now, we are trying to minimize 



E = J2 (* t + Ci • 6b + di ■ Suf , 



t=i 
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subject to the condition b • 8b = 0. This constrained minimization problem 
can be transformed into an equivalent unconstrained form by introduction 
of a Lagrangian multiplier. Instead of minimizing E itself, we then have to 
minimize: 

E' = E + A(b • (5b). 

Differentiating E' with respect to 8b, and setting the result equal to zero 
yields, 

n 

J^ (ti + Ci ■ 8b + dj • 8u) Cj + Ab = 0. 

t=i 
Before we can proceed, we need to eliminate the unknown Lagrangian multi- 
plier A. The dot-product of the expression with b leads to 

n 

Y^ (U + Ci • 8b + d; • 6ui) a • b + A(b • b) = 0, 

t=i 
which, since b • b = 1, gives us a value for A that we can use to compute the 
term 

n 

Ab = — b 2_. (ti + c,- • <!>b + dj • 8u) c* • b, 

i=l 

or 



Ab = -(bb r ) ]T (U + ci ■ 8b + di ■ 8u) a. 
Finally, substituting for Ab in the equation above we obtain 

n 

B^2(U + Ci 8b + dj 6u) Ci = 0, 

i-l 

where 

B = (I- bb T ) 

is the projection operator that removes components of vectors parallel to the 
baseline, and i" is just the 3x3 identity matrix. We can conclude that the 
equation above relates quantities in a plane perpendicular to the baseline. 

Finally, if we differentiate E' with respect to 8w and set this result also 
equal to zero, we obtain 

n 

J^ (U + Ci ■ 8b + di ■ 6u) di = 0. 
i=i 
Together, the two vector equations constitute six linear scalar equations in 
the six unknown components of 8b and 8u>. We can rewrite them in the more 
compact form: 

BC 8b + BF 8u = -Bc 
F T 8b + D 8u> = -d 
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or 



where 



(BC BF\ ( 8h\ _ _ ( Bc\ 

n n n 

C = J^c,-cf, F = Y J c * d L and £> = X] d ' d ^' 

t=l «=1 i=l 



W 



hile 



c = YijCj and d — y tjdj. 

i=l i-1 

The matrix B is singular, as can be seen by noting that b is an eigenvector with 
zero eigenvalue. Thus the first three equations above are not independent. 
One of them will have to be removed. Fortunately, we also still have to 
incorporate the condition that 6b be perpendicular to b. We can do this by 
replacing one of the first three equations with the linear equation b • 6b = 0. 
For the best numerical accuracy one should eliminate the equation with the 
smallest coefficients. 

The above gives us a way of finding small changes in the baseline and 
rotation that reduce the overall error sum. This method can be applied itera- 
tively to locate a minimum. Numerical experiments confirm that it converges 
rapidly when a good initial guess is available. Incremental adjustments can- 
not be computed if the six-by-six coefficient matrix becomes singular. This 
will happen when there are fewer than five pairs of rays, and for certain rare 
configurations of points in the scene (see the discussion of critical surfaces 
later on). 

7. Adjusting the Baseline and the Rotation 

The iterative adjustment of the baseline is straightforward: 

b n+l =h n +6h n^ 

where b n is the baseline estimate at the beginning of the n-th iteration, while 
6b n is the adjustment computed during the n-th iteration, as discussed in 
the previous section. If <!>b n is not infinitesimal, the result will not be a unit 
vector. We can, and should, normalize the result by dividing by its magnitude. 

7.1. Adjustment of Rotation using Unit Quaternions 

Adjusting the rotation is a little harder. Rotations are conveniently repre- 
sented by unit quaternions [Stuelpnagle 1964, Salamin 1979, Taylor 1982, 
Horn 1986, Horn 1987a]. The groundwork for the application of the unit 



Relative Orientation 13 

quaternion notation in photogrammetry was laid by Thompson [1959a] and 
Schut [1958-59]. A positive rotation about the axis u; through an angle 6 is 
represented by the unit quaternion 

q = cos(0/2) + sin(0/2) w, 

where ui is assumed to be a unit vector. Composition of rotations corresponds 
to multiplication of the corresponding unit quaternions. The rotated version 
of a vector r is computed using 

o/ o o o * 

r =qrq , 

where q* is the conjugate of the quaternion q , that is, the quaternion obtained 
by changing the sign of the vector part. Here, f is a purely imaginary quater- 
nion with vector part r, while ?' is a purely imaginary quaternion with vector 
part r'. 

The infinitesimal rotation 8u corresponds to the quaternion 

8uj = 1 + 8u). 

We can adjust the rotation q by premultiplying with 8u, that is, 

q n+1 = 8u n q n . 

If 8u3 n is not infinitesimal, 8Cj h will not be a unit quaternion, and so the result 
of the adjustment will not be a unit quaternion either. This undesirable state 
of affairs can be avoided by using either of the two unit quaternions 

Su = y/l- HHI 2 + to, 
or 

8u> = (l+8u>)/Jl + \\8u>\\ 2 . 

Alternatively, one can simply normalize the product by dividing by its mag- 
nitude. 

7.2. Adjustment of Rotation using Orthonormal Matrices 

The adjustment of rotation is a little trickier if orthonormal matrices are used 
to represent rotations. We can write the relationship 

r' = r + (8u x r), 

in the form 

r' = r + W r, 

where the skew-symmetric matrix W is defined by 

— 8(jJ Z 8(jJy 

W = | 8u z -8u) x | , 

— 8u>y 8li) x 
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in terms of the components of rotation vector 8u = (8(v x ,8u} y ,8u} z ) T . Conse- 
quently we may write r' = Qr, where Q — I + W, or 

(1 —8u z Suiy 
8u z 1 —8u x 

—8u) y 8uj x 1 

One could then attempt to adjust the rotation by multiplication of the ma- 
trices Q and R as folows: 

R n+1 = Qn R n 

The problem is that Q is not orthonormal unless 8w is infinitesimal. In prac- 
tice this means that the rotation matrix will depart more and more from 
orthonormality as more and more iterative adjustments are made. It is pos- 
sible to re-normalize this matrix by finding the nearest orthonormal matrix, 
but this is complicated, involving the determination of the square-root of a 
symmetric matrix [Horn, Hilden & Negahdaripour 1988]. This is one place 
where the unit quaternion representation has a distinct advantage: it is trivial 
to find the nearest unit quaternion to a quaternion that does not have unit 
magnitude. 

To avoid this problem, we should really start with an orthonormal ma- 
trix to represent the incremental rotation. We can use either of the unit 
quaternions 

8u = y/l - \\8l>\\ 2 + Sta, 
or 

6u> = (l + 6u,)/y/l + \\6u\\ 2 , 
to construct the corresponding orthonormal matrix 

(00 + 9* - 1 2 y ~ <& 2(<lx<ly ~ <k>qz) 2(q x q z + q q y ) 

2(q y qx + qoqz) ql-ql + q\- ql i{q y q z - q q x ) 
2(qzqx - qoq y ) 2(q z q y + q q x ) ql - q\ - q 2 y + q\ 

where qo is the scalar part of the quaternion 8u>, while q x , q y , q z are the com- 
ponents of the vector part. Then the adjustment of rotation is accomplished 
using 

R n+1 = Qn R n 

Note, however, that the resulting matrices will still tend to depart somewhat 
from orthonormality due to numerical inaccuracies. This may be a problem 
if many iterations are required. 

8. Inherent Ambiguities 

The iterative adjustment described above may arrive at a number of appar- 
ently different solutions. Some of these are just different representations of 
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the same solution, while others are related to the correct solution by a simple 
transformation. First of all, note that — q represents the same rotation as q, 
since 

-q)r(-q ) = qrq . 
That is, antipodal points on the unit sphere in four dimensions represent the 
same rotation. We can prevent any confusion by ensuring that the first non- 
zero component of the resulting unit quaternion is positive, or that the largest 
component is positive. 

Next, note that the triple product, [b r{ r>], changes sign, but not mag- 
nitude, when we replace b with — b. Thus the two possible senses of the 
baseline yield the same sum of squares of errors. However, changing the sign 
of b does change the signs of both a and /?. All scene points imaged are in 
front of the camera, so the distances should all be positive. In the presence 
of noise, it is possible that some of the distances turn out to be negative, but 
with reasonable data almost all of them should be positive. This allows us to 
easily pick the correct sense for the baseline. 

Not so obvious is another possibility: Suppose we turn all of the left mea- 
surements through 7r radians about the baseline, in addition to the rotation 
already determined. That is, replace r{ with 

r',' = 2(b • r{)b - r}. 
This follows from Rodrigues' formula for the rotation of a vector r: 

cos 6 r + sin 8 (u> x r) + (1 — cos 9)(u3 ■ r) u>, 
with 6 — 7r and u> = b. Then the triple product [b rj r>] turns into 

2(b • r',)[b b r r ] - [b r{ r r ] = -[b rj r r ]. 
This, once again, reverses the sign of the error term, but not its magnitude. 
Thus the sum of squares of errors is unaltered. The signs of a and /? are 
affected, however, although this time not in as simple a way as when the 
sense of the baseline was reversed. 

If [b rj r r ] = 0, we find that exactly one of a and /3 changes sign. This 
can be shown as follows: The triple product will be zero when the left and 
right rays are coplanar with the baseline. In this case we have 7 = 0, and so 

ar't = b + /?r r , 
Taking the cross-product with b we obtain 

a(r',xb) = /?(r r xb), 
If we now replace rj by r" = 2(b • rj)b — rj, we have for the new distances a' 
and /3' along the rays: 

-c/(r; xb) = /?'(r r xb), 
We conclude that the product a' /?' has sign opposite to that of the product 
a/3. So if a and /? are both positive, one of a' or j3' must be negative. 



16 

In the presence of measurement error the triple product will not be ex- 
actly equal to zero. If the rays are nearly coplanar with the baseline, however, 
we find that one of a and /3 almost always changes sign. With very poor data, 
it is possible that both change sign. (Even with totally random ray directions, 
however, this only happens 27.3% of the time, as determined by Monte Carlo 
methods). In any case, we can reject a solution in which roughly half the 
distances are negative. Moreover, we can find the correct solution directly by 
introducing an additional rotation of 7r radians about the baseline. 

9. Remaining Ambiguity 

If we take care of the three apparent two-way ambiguities discussed in the 
previous section, we find that in practice a unique solution is found, provided 
that a sufficiently large number of ray pairs are available. That is, the method 
converges to the unique global minimum from every possibly starting point 
in parameter space. 

Local minima in the sum of squares of errors appear when only a few 
more than the minimum of five ray pairs are available (as is common in prac- 
tice). This means that one has to repeat the iteration with different starting 
values for the rotation in order to locate the global minimum. A starting 
value for the baseline can be found in each case using the closed-form method 
described in section 5. To search the parameter space effectively, one needs 
a way of efficiently sampling the space of rotations. The space of rotations 
is isomorphic to the unit sphere in four dimensions, with antipodal points 
identified. The rotation groups of the regular polyhedra provide convenient 
means of uniformly sampling the space of rotations. The group of rotations of 
the tetrahedron has 12 elements, that of the hexahedron and the octahedron 
has 24, and that of the icosahedron and the dodecahedron has 60 (representa- 
tions of these groups are given in Appendix A for convenience). One can use 
these as starting values for the rotation. Alternatively, one can just generate 
a number of randomly placed points on the unit sphere in four dimensions as 
starting values for the rotation. 

When there are only five pairs of rays, the situation is different again. In 
this case, we have five non-linear equations in five unknowns and so in general 
expect to find a finite number of exact solutions. That is, it is possible to 
find baselines and rotations that satisfy the coplanarity conditions exactly 
and reduce the sum of squares of errors to zero. If we ignore the three types 
of ambiguities discussed above (incuding the rotation by tt radians about the 
baseline), then there are generally four distinct sets of baselines and rotations 
that yield an exact solution. Typically only one of these yields positive signs 
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for all of the distances. These are empirical observations; I have not been able 
to prove that there are generally four solutions that satisfy the coplanarity 
conditions. 

10. Summary of the Algorithm 

Consider first the case where we have an initial guess for the rotation. We 
start by finding the best-fit baseline direction using the closed-form method 
described in section 5. We determine the correct sense of the baseline by 
choosing the one that makes most of the signs of the distances positive. Then 
we proceed as follows: 

• For each pair of corresponding image points, we compute rj t , the left 
ray direction r/^ rotated into the right camera coordinate system, 
using the present guess for the rotation. 

• We then compute the cross-product c, = r{ i x r r> i, the double cross- 
product di = r'f t x (r rj i x b) and the triple-product t{ = [b r\ i r ri i\. 

• We accumulate the dyadic products c,cf, C{df and didf, as well 
as the vectors t,c,- and t{di. The totals of these quantities over all 
image point pairs give us the matrices C, F, D and the vectors c 
and d. 

• We can now solve for the increment in the baseline Sh and the in- 
crement in the rotation 8u> using the method derived in section 6. 

• We adjust the baseline and the rotation using the methods discussed 
in section 7, and recompute the sum of the squares of the error terms. 

The new orientation parameters are then used in the next iteration of the 
above sequence of steps. As is the case with most iterative procedures, it is 
hard to decide when to stop. Typically, the total error becomes small after a 
few iterations and no longer decreases at each step, because of limited accuracy 
in the arithmetic operations. One can stop the iteration the first time the error 
increases. Alternatively, one can stop after either a fixed number of iterations 
or after the error becomes less then some predetermined threshold. 

When the decision has been made to stop the iteration, a check of the 
signs of the distances along the rays is in order. If most of them are negative, 
the baseline direction should be reversed. If neither sense of the baseline 
direction yields mostly positive distances, one needs to consider the possibility 
of a rotation through tt radians about the baseline b. If this also yields mixed 
signs, one is dealing with a local extremum of the error function; something 
that will only happen if the initial guess is in fact not adequate. 

If an initial guess is not available, one proceeds as follows: 
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• For each rotation in the chosen group of rotations, perform the above 
iteration to obtain a candidate baseline and rotation. 

• Choose the solution that yields the smallest total error. 

When there are many pairs of rays, the iterative algorithm will converge to 
the global minimum error solution from any initial guess for the rotation. 
There is no need to sample the space of rotations in this case. Also, instead 
of sampling the space of rotations in a systematic way using a finite group 
of rotations, one can generate points randomly distributed on the surface of 
the unit sphere in four-dimensional space. This provides a simpler means of 
generating initial guesses, although more initial guesses may have to be tried 
than when a systematic procedure is used to sample the space of rotations. 
The method given minimizes the sum of the squares of the triple products 

[b rj r r ]. 

If desired, one can modify it to use some multiple of the triple product as an 
error term by weighting the contribution to the overall sum. This can lead to a 
considerably more complex optimization problem if the weights depend on the 
unknown baseline and the unknown rotation. This happens, for example, if 
we try to minimize the sum of squares of the sines of the angles corresponding 
to vertical disparity: 

«=JM 



lib x ri|| ||b xr r ||- 

If we assume that the weighting factors vary slowly during the iterative pro- 
cess, however, we can to use the current estimates of the baseline and rotation 
in computing the weighting factors, and not take into account the small vari- 
ations in the error sum due to changes in the weighting factors. That is, when 
taking derivatives, the weighting factors are considered constant. This is a 
good approximation when the parameters vary slowly, as they will when one 
is close to an extremum. 

11. Critical Surfaces 

In certain rare cases, relative orientation cannot be recovered fully, even when 
there are five or more pairs of rays. Normally, each error term varies linearly 
with distance in parameter space from the location of an extremum, and so the 
sum of squares of errors varies quadratically. There are situations, however, 
where the error terms to not vary linearly with distance, but quadratically or 
higher order, in certain special directions in parameter space. In this case, the 
sum of squares of errors does not vary quadratically with distance from the 
extremum, but as a function of the fourth or higher power of this distance. 
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This makes it very difficult to accurately locate the extremum. In an extreme 
situation, the total error may not vary at all along some curve in parameter 
space. In this case, the total error is unaffected by a small change in the 
rotation, as long as this change is accompanied by a corresponding small 
change in the baseline. There is no localized extremum and consequently 
the solution for relative orientation is not unique. It turns out that this 
problem arises only when the observed scene points lie on certain surfaces 
called Gefdhrliche Fldchen or critical surfaces [Brandenberger 1947, Hofmann 
1949, Zeller 1952, Schwidefsky 1973]. We show next that only points on 
certain hyperboloids of one sheet and their degenerate forms can lead to this 
kind of problem (see also [Horn 1987b]). 

We could try to find a direction of movement in parameter space (£b, 6u>) 
that leaves the total error unaffected when given a particular surface. Instead, 
we will take the critical direction of motion in the parameter space as given, 
and try to find a surface for which the total error is unchanged. 

Let R be a point on the surface, measured in the right camera coordinate 
system. Then 

/3 r r = R and a rj = b + R, 
for some positive a and /?. In the absence of measurement errors, 

[brJr r ] = i[b(b + R)R]=0. 

We noted earlier that when we change the baseline and the rotation slightly, 
the triple product [b rj r r ] becomes 

[(b + flj)(r j , + &)xrl)r r ] ) 

or, if we ignore higher-order terms, 

[b r{ r r ] + (r{ x r r ) • 6b + (r r x b) • (Su x rj). 

The problem we are focusing on here arises when this error term is unchanged 
for small movement in some direction in the parameter space. That is when 

(rj x r r ) • Sb + (r r x b) • (<W x r{) = 0, 

for some 6b and 6ui. Introducing the coordinates of the imaged points we 
obtain: 

-^ (((b + R) x R) • <5b + (R x b) • (Su x (b + R))) = 0, 

or 

(R x b) • (Su x R) + (R x b) • (6e> x b) + [b R Sb] = 0. 

If we expand the first of the dot-products of the cross-products, we can write 
this equation in the form 

(R • b)(6u • R) - (b • Su)(R R) + L R = 0, 
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where 



L = fxb, while I = b x 6<v + Sb. 



The expression on the left-hand side contains a part that is quadratic in R 
and a part that is linear. The expression is clearly quadratic in X, Y, and 
Z, the components of the vector R = (X,Y,Z) T . Thus a surface leading to 
the kind of problem described above must be a quadric surface [Korn & Korn 
1968]. 

Note that there is no constant term in the equation of the surface, so 
R = satisfies the equation. This means that the surface passes through the 
right projection center. It is easy to verify that R = — b satisfies the equation 
also, which means that the surface passes through the left projection center 
as well. In fact, the whole baseline (and its extensions), R = fcb, lies in the 
surface. This means that we must be dealing with a ruled quadric surface. It 
can consequently not be an ellipsoid or hyperboloid of two sheets, or one of 
their degenerate forms. The surface must be a hyperboloid of one sheet, or 
one of its degenerate forms. Additional information about the properties of 
these surfaces is given in Appendix B, while the degenerate forms are explored 
in Appendix C. 

It should be apparent that this kind of ambiguity is quite rare. This is 
nevertheless an issue of practical importance, however, since the accuracy of 
the solution is reduced if the points lie near some critical surface. A textbook 
case of this occurs in aerial photography of a roughly U-shaped valley taken 
along a flight line parallel to the axis of the valley from a height above the 
valley floor approximately equal to the width of the valley. In this case the 
baseline lies on a circular cylinder that also lies close to the surface on which 
the imaged points lie. This means that it is close to one of the degenerate 
forms of the hyperboloid of one sheet (see Appendix C). 

Note that hyperboloids of one sheet and their degenerate forms are ex- 
actly the surfaces that lead to ambiguity in the case of motion vision. The 
coordinate systems and symbols have been chosen here to make the correspon- 
dence between the two problems more obvious. The relationship between the 
two situations is nevertheless not quite as transparent as I had thought earlier 
[Horn 1987b]. In the case of the ambiguity of the motion field, for example, we 
are dealing with a two-way ambiguity arising from infinitesimal displacements. 
Here we are dealing with an infinite number of solutions arising from images 
taken with cameras that have large differences in position and orientation. 
Also note that the symbol Siv stands for a small change in a finite rotation 
here, while it refers to a difference in instantaneous rotational velocities in the 
motion vision case. 
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12. Conclusions 

Methods for recovering the relative orientation of two cameras are of impor- 
tance in both binocular stereo and motion vision. A new iterative method for 
finding the relative orientation has been described here. It can be used even 
when there is no initial guess available for the rotation or the baseline. The 
new method does not use Euler angles to represent the orientation. 

When there are many pairs of corresponding image points, the iterative 
method finds the global minimum from any starting point in parameter space. 
Local minima in the sum of squares of errors occur, however, when there 
are relatively few pairs of corresponding image points available. Method for 
efficiently locating the global minimum in this case were discussed. When only 
five pairs of corresponding image points are given, several exact solutions of 
the coplanarity equations can be found. Typically only one of these yields 
positive distances to the points in the scene, however. This allows one to pick 
the correct solution even when there is no initial guess available. 

All of these methods fail in the rare case that the scene points lie on a 
critical surface. 
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Appendix A — Rotation Groups of Regular Polyhedra 

Each of the rotation groups of the regular polyhedra can be generated from 
two judiciously chosen elements. For convenience, however, an explicit rep- 
resentation of all of the elements of each of the groups is given here. The 
number of different component values occuring in the unit quaternions rep- 
resenting the rotations can be kept low by careful choice of the alignment of 
the polyhedron with the coordinate axes. The attitudes of the polyhedra here 
were selected to minimize the number of different numerical values that occur 
in the components of the unit quaternions. A different representation of the 
group is obtained if the vector parts of each of the unit quaternions is rotated 
in the same way. This just corresponds to the rotation group of the polyhe- 
dron in a different attitude with respect to the underlying coordinate system. 
This observation leads to a convenient way of generating finer systematic sam- 
pling patterns of the space of rotations than the ones provided directly by the 
rotation group of a regular polyhedron in a particular alignment with the 
coordinate axes (see also [Brou 1983]). 

The components of the unit quaternions here may take on the values 
and 1, as well as the following: 

y/B-1 . 1 1 , , y/E+l 

a=—^—, 6=-, c=—, and d=—^—. 

Here are the unit quaternions for the twelve elements of the rotation 
group of the tetrahedron: 

(1, 0, 0, 0) (0, 1, 0, 0) (0, 0, 1, 0) (0, 0, 0, 1) 
(6, 6, 6, 6) (b, b, 6,-6) (6, 6,-6, 6) (6, 6,-6,-6) 
(6,-6, 6, 6) (6,-6, 6,-6) (6,-6,-6, 6) (6,-6,-6,-6) 

Here are the unit quaternions for the twenty-four elements of the rotation 
group of the octahedron and the hexahedron (cube): 

(1, 0, 0, 0) (0, 1, 0, 0) (0, 0, 1, 0) (0, 0, 0, 1) 

(0, 0, c, c) (0, 0, c,-c) (0, c, 0, c) (0, c, 0,-c) 

(0, c, c, 0) (0, c,-c, 0) (c, 0, 0, c) (c, 0, 0,-c) 

(c, 0, c, 0) (c, 0,-c, 0) (c, c, 0, 0) (c,-c, 0, 0) 

(6, 6, 6, 6) (6, 6, 6,-6) (6, 6,-6, 6) (6, 6,-6,-6) 

(6,-6, 6, 6) (6,-6, 6,-6) (6,-6,-6, 6) (6,-6,-6,-6) 

Here are the unit quaternions for the sixty elements of the rotation group of 
the icosahedon and the dodecahedron: 
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(d, 6, -a, 0) 


(d, -6, a, 0) 


(d,-6,-a, 0) 


b) 


(6, 6, 6,-6) 


(6, 6,-6, 6) 


(ft, ft, -ft, -b) 


b) 


(6,-6, 6,-6) 


(6,-6,-6, 6) 


(6,-6,-6,-6) 



Remember that changing the signs of all the components of a unit quaternion 
does not change the rotation that it represents. 



Appendix B — Some Properties of Critical Surfaces 

In this appendix we develop some more of the properties of the critical sur- 
faces. The equation of a critical surface can be written in the form 

(R x b) • (6w x R) + L • R = 0, 

or 

(R • b)(Sio • R) - (b • 6w)(R R) + L R = 0, 
where 

L = £ x b, while £ = b x Sw + Sb. 
It is helpful to first establish some simple relationships between the quantities 
appearing in the formula above. We start with the observations that £• b = 0, 
that £ ■ 8u = 6b ■ Su, and £ x Sb = —{Sb • 8u))b. 
We can also expand L to yield, 

L = 6u-(b- 8u)b + Sb x b. 
It follows that L • b = 0, that L • Sb = 8w ■ Sb, and 

L x b = -(b x 8u + Sb) = -£, 
Lx8u> = (8b- 8u)b - (b • 6u)L 

We have already established that R = k b is an equation for one of the rulings 
passing through the origin. A hyperboloid of one sheet has two intersecting 
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families of rulings, so there should be a second ruling passing through the 
origin. Consider the vector S denned by 

S = (L x Su) x L, 
which can be written in the form 

S = (L • L)Su - (L • Su)L, 
or 

S = (6u>- Sb)£ + (b • Su)(£ • £)b, 
so that S • b = (b • Su)(£ ■ £) and S • Su = (Su ■ 8b) 2 + (b • Su) 2 (£ ■ £). 
If we substitute R = k S into the formula 

(R x b) • (Su x R) + L • R, 
we obtain zero, since L • S = and 

S x b = (Su • Sb)L 
is orthogonal to 

S x Su = — (L • Su) L x Su. 
We conclude that R = A;S is an equation for the other ruling that passes 
through the right projection center. 

There are two families of parallel planes that cut an ellipsoid in circular 
cross-sections [Hilbert & Cohn-Vossen 1953]. Similarly, there are two families 
of parallel planes that cut a hyperboloid of one sheet in circular cross-sections. 
One of these families consists of planes perpendicular to the baseline, that is, 
with common normal b. We can see this by substituting R • b = k in the 
equation of the critical surface. We obtain 

k(Su • R) - (b • 8u)(R ■ R) + L • R = 0, 
or 

(b • 8u)(TL ■R)-(k8u + L)-R = 0. 
This is the equation of a sphere, since the only second-order term in R is a 
multiple of 

R R = X 2 +Y 2 + Z 2 . 
We can conclude that the intersection of the critical surface and the plane is 
also the intersection of this sphere and the plane, and so must be a circle. 
The same applies to the intersection of the critical surface and the family of 
planes with common normal Su, since we get 

(b • 6u)(R ■ R) - (Jfe b + L) • R = 0, 
when we substitute R • Su = k into the equation of the critical surface. 

The equation of the critical surface is given in the implicit form /(R) = 
0. The equation of a tangent plane to such a surface can be obtained by 
differentiating with respect to R: 

N = (Rx^u)xb + (Rxb)x^+L 
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or 

N = (R • b)Su + (R • 6u})b - 2(b • 8u>)R + L 
The tangent plane at the origin has normal L. This tangent plane contains 
the baseline (since L • b = 0), as well as the other ruling passing through the 
origin (since L • S = 0). Note that the normal to the tangent plane is not 
constant along either of these rulings, as they would be if we were dealing 
with a developable surface. 

In the above we have not considered a large number of degenerate situ- 
ations that can occur. The reader is referred to Appendix C for a detailed 
analysis of these. 

Appendix C — Degenerate Critical Surfaces 

There are a number of special alignments of the infinitesimal change in the 
rotation, 8u>, with the baseline, b, and the infinitesimal change in the baseline 
8b that lead to degenerate forms of the hyperboloid of one sheet. 

One of the rulings passing through the origin is given by R = k b, while 
the other is given by R = k S. If these two rulings become parallel, we are 
dealing with a degenerate form that has only one set of rulings, that is a 
conical surface. Now 

S = (6u- 6h)£ + (b • 6u)(l ■ £)b, 

is parallel to b only when (8u> • 8b) = 0, since I is perpendicular to b. In this 
case 

8b ■ 8u = and 8b b = 0, 
so 8b = A;(b x 8ui) for some constant k. Consequently I = (k + l)(b x 8u), 
The vertex of the conical surface must lie on the baseline since the baseline 
is a ruling, and every ruling passes through the vertex. It can be shown that 
the vertex actually lies atR=— (k + l)b. 

We also know that cross-sections in planes perpendicular to the baseline 
are circles. This tells us that we are dealing with elliptical cones. Right 
circular cones cannot be critical surfaces. It can be shown that the main axis 
of the elliptical cone lies in the direction b + 8ut. 

A special case of the special case above occurs when 

||b x 8u>\\ = 0, 

that is 8lo || b. Here 8u = kb for some constant k and so I = 8b and 
Li — 8b x b. The equation of the surface becomes 

fc(R • b) 2 - Jb(b • b)(R • R) + L • R = 0, 
or 

k ||R x b|| 2 + (8b x b) • R = 0. 
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This is the equation of a circular cylinder with axis parallel to the baseline. 
In essence, the vertex of the cone has receded to infinity along the baseline. 

Another special case arises when the radius of the circular cross-sections 
with planes perpendicular to the baseline becomes infinite. In this case we 
obtain straight lines, and hence rulings, in these planes. The hyperbolic 
paraboloid is the degenerate form that has the property that each of its two 
sets of rulings can be obtained by cutting the surface with a set of parallel 
planes [Hilbert & Cohn-Vossen 1953]. This happens when 8u is perpendicular 
to b, that is, b • 8u = 0. The equation of the surface in this case simplifies to 

(R • b)(8ui R) + L R = 0. 

The intersection of this surface with any plane perpendicular to the baseline 
is a straight line. We can show this by substituting 

Rb = k, 

into the equation of the surface. We obtain 

(Jfc*w + L)-R = 0, 
that is, the equation of another plane. Now the intersection of two planes is 
a straight line. So we may conclude that the intersection of the surface and 
the original plane is a straight line. We can show in the same way that the 
intersection of the surface with any plane perpendicular to 8u is a straight 
line by substituting 

R • 8uj = k 
into the equation of the surface. It can be shown that the saddle point of the 
hyperbolic paraboloid surface lies on the baseline. 

A special case of particular interest arises when 8u is perpendicular to 
both b and 8b, that is, 

b • 8u — and 8b • 8u> = 0, 

and so 8u = k(8b x b), for some constant k. Then £ — (k + l)8b and 
L = (k + l)(£b x b). The equation of the surface becomes 

fc(R • b)((8b xb)-R)+(H l)((8b x b) • R) = 0, 
or just 

(k(K -b) + (k + 1)) ((8b x b) • R) = 0. 
so either 

(8b x b) • R = or Jb(R • b) + (Jb + 1) = 0. 
The first of these is the equation of a plane containing the baseline b and 
the vector 8b. The second is the equation of a plane perpendicular to the 
baseline. So the solution degenerates in this case into a surface consisting of 
two intersecting planes. One of these planes appears only as a line in each of 
the two images, since it passes through both projection centers, and so does 
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not really contribute to the image. It is fortunate that planes can only be 
degenerate surfaces if they are perpendicular to the baseline, since surfaces 
that are almost planar occur frequently in aerial photography. 1 
To summarize then, we have the following degenerate cases: 

• elliptical cones when 8u) J_ £b, 

• circular cylinders when 8u || b, 

• hyperbolic parabloids when 8(v ± b, 

• intersecting planes when 8u ± 8b and 8u ± b. 

For further details, and a proof that not all hyperboloids of one sheet passing 
through the origin lead to critical surfaces, see [Horn 1987b]. 



The baseline was nearly perpendicular to the surface in the sequence of pho- 
tographs obtained by the Ranger spacecraft as it crashed into the lunar surface. 
This made photogrammetric analysis difficult. 
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